Residual Coulomb interaction fluctuations in chaotic systems: the boundary, random plane waves, 

and semiclassical theory 
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New fluctuation properties arise in problems where both spatial integration and energy summation are nec- 
essary ingredients. The quintessential example is given by the short-range approximation to the first order 
ground state contribution of the residual Coulomb interaction. The dominant features come from the region 
near the boundary where there is an interplay between Friedel oscillations and fluctuations in the eigenstates. 
Quite naturally, the fluctuation scale is significantly enhanced for Neumann boundary conditions as compared 
to Dirichlet. Elements missing from random plane wave modeling of chaotic eigenstates lead surprisingly to 
significant errors, which can be corrected within a purely semiclassical approach. 
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The characterization of quantum systems with any of a 
variety of underlying classical dynamics, ranging from dif- 
fusive to chaotic to regular, has often demonstrated that the 
study of their statistical properties is of primary importance. 
Spectral fluctuations are a principle example as they gave the 
first support to one of the main results linking classical chaos 
and random matrix theory Hf], the Bohigas-Giannoni-Schmit 
conjecture fH Needless to say, the statistical properties 
of eigenfunctions are also a subject of paramount interest 

lilaiiEaOiii. 



For chaotic systems, a widely accepted starting point for 
the treatment of eigenfunction fluctuations locally, such as 
the amplitude distribution or the two-point correlation func- 
tion c(|r — r'|) = (ip(r)ip(r')) °f a gi ven eigenfunction t/j, is 
a model ing in terms of a random superposition of plane waves 
(RPW) lUBI. For two-degree-of-freedom systems, c(r) can 
be understood as being given approximately by a Bessel func- 
tion. For distances |r — r'| short compared to the system size, 
and in the absence of effects related to classical dynamics 
[Till , this is roughly observed in numerical J3, FA [l5 1 
and experimental Hill studies. Our interest here is in statisti- 
cal properties of eigenfunctions going beyond local quantities 
such as c(r). 

One motivation for the introduction of these new statisti- 
cal measures is to study the interplay between interferences 
and interactions in mesoscopic systems. For typical elec- 
tronic densities, the screening length is close to the Fermi 
wavelength A F , and the screened Coulomb interaction can be 
approximated by the short range expression V sc (r — r') = 
(2v)~ 1 Fq5(t — r') with 2v the mean local density of states, 
including spin degeneracy, (v = m/2Tih 2 for d = 2) and F£ 
the dimensionless Fermi liquid parameter lfl7ll . a constant of 
order one. 



To this level of approximation, the first order ground state 
energy contribution of the residual interactions can be ex- 
pressed in the form 6E RI = (2v)- 1 Fg J dr N^(r)Ni(r), 
with N a the unperturbed ground state density of particles with 
spin a. From this expression, it is seen that the increase of in- 
teraction energy associated with the addition of an extra elec- 
tron is related to 



Si = A 



fdr\%(r)\ 2 N(r;E+), (1) 



with N(r; E) = J2Zi l*i( r )l e ( E ~ E i) and the under- 
standing that Ei < Ef < Ei + i. Our goal in this letter is 
to study the fluctuation properties of the Si, concentrating on 
the case of two dimensional billiards with either Dirichlet or 
Neumann boundary conditions. For them, A is the billiard 
area. 

The dominant contributions to Si (and its fluctuations) orig- 
inate from the Friedel oscillations of the density of particles 
near the boundary. For billiards systems, they can be ex- 



pressed as iV Friedel (r; E) 



N W (E) 
A 



1± 



[18], with 



J x (2kx) 
kx 

x the distance from the boundary, the + and — sign corre- 
sponding respectively to Neumann and Dirichlet boundary 
conditions, and Nw(E) refers to leading term of the Weyl 
formula, Nw(E) = vAE. To leading order we can therefore 
use the approximation 



Si = i ± i / dr 



J\(2k F x) 



*i(r)| 5 



(2) 
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To proceed, a description of the fluctuations of |5'i(r)| is 
also required. These are obtained ahead using a semiclassi- 
cal approach closely related to the Gutzwiller trace formula. 
However, it is useful first to consider the oft-employed RPW 
description, which very interestingly turns out to lack a cou- 
ple of crucial ingredients. Nevertheless, it sheds light on the 
mechanism governing the fluctuations under study. 

Within RPW |4l |5[] eigenstates are represented, in the ab- 
sence of any symmetry, by a random superposition of plane 
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waves J^i a i exp(ik/ • r) with wave-vectors of fixed modulus 
|k/| = k F distributed isotropically. Time reversal invariance 
introduces a correlation between time reversed plane waves 
such that the eigenfunctions are real. Similarly, the presence 
of a planar boundary imposes a constraint between the coeffi- 
cients of plane waves related by a sign change of the normal 
component of the wave-vector k; \19, 2(J • Near a boundary, 
and using a system of coordinates r = x + y with (x,y) the 
vectors respectively parallel and perpendicular to the bound- 
ary, eigenfunctions are mimicked statistically by a superposi- 
tion, 
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— — V"a;cs(k ; -x)cos(k; ■ f + <fil) (3) 
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where cs(-) = f sin(-) for Dirichlet and cos(-) for Neumann 
boundary conditions. The phase angle <pi, the orientation of 
the wave vector k;, and the real amplitude a; with ifliai') — 
5u'(j 2 are all chosen randomly. Normalization of the wave- 
functions fixes the variance a 2 through the relation 



dr (l^(r)l' 



Aa 2 
4JV, 



1 ± 



C 



2k F A 



(4) 



where C is the perimeter. 

The variance is a natural measure of the fluctuations. To 
leading order 



Var (S t ) 



dridr- 



Ji(2k F Xi) J\(2k F X2) 



A 2 J 1 k F xi k F X2 

x [(|^(r 1 )| 2 |* i (r 2 )| 2 )-(|^(r 1 )| 2 )(|* l (r 2 )| 2 )] . 

The fluctuations are thus given by pair-wise correlating the 
random plane wave coefficients a n such that (a/a// a m a m i) = 
(o/a//) (a m a m >) + (aja m ) {a v a m ,) + (aia m i){ai>a m ) = 
cr 4 (Su'Smm' + Simh'm' + <W<W)- One obtains in this way 



Var [Si 



8i 2 1 



AT eff 
eff l,m=l 



dridr 2 



Ji{2k F x\) J\(2k F X2) 



k F x 2 



cs(k; • xi)cs(k; • x 2 )cs(k m • xi)cs(k m • x 2 ) 
cos [ki ■ (yi - y 2 )] cos [k m • (yi - y 2 )] . 

Performing the integral gives 

VarN = ^<A 2 (#)) e 



(5) 



(6) 



(lef 

for 



where we have introduced the function X(9) 
[l ± | sin(0)|] /| cos(0)| and the expression 

(\ 2 (0)) e = d(sin0)\ 2 (8) is given ahead in Eq. O- 

In Fig. [1] the variance of Si as a function of i is repre- 
sented for a chaotic system, the cardioid billiard with Dirichlet 
boundary conditions 12 ll 12211 . In this case, quite surprisingly 
given the history of modeling chaotic eigenstates within the 
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FIG. 1 : Variance of Si for the cardioid billiard with Dirichlet bound- 
ary conditions. The solid line is the RPW expression, Eq. (6), the 
dashed line is the semiclassical prediction, Eq. d 1 4b . and the discrete 
points are for the cardioid billiard (shown in inset). 



RPW framework, its predictions significantly overestimate the 
fluctuations. In fact, two important elements are missing from 
this approach, both of which are addressed properly within a 
purely semiclassical approach. 

One difficulty immediately encountered with a semiclas- 
sical treatment of the Si is that their computation im- 
plies addressing the fluctuations of individual wave-functions, 
whereas the semiclassical approximations valid for chaotic 
systems of use here converge only for (locally) smoothed 
quantities. This difficulty may be overcome by following the 
spirit of Bogomolny's calculation 02311 . and introducing a local 
energy averaging, 



(S) 



del' 



AN 



— V 

AN ^ 

E -AE <El< E+^f 



(7) 



with AN = N(E+2f)-N(E-2f). This generates the 
relation 



Var[(5) A Ar] d = ((S)an-S)' 



AN 



Var^; 



1 



AN J 



Covar [SiS., 



i+j) 



(8) 



assuming that there is translational invariance (smooth, uni- 
form behavior) locally in indices (i, j). See Il24ll for further 
discussion for the relevance of this assumption. Computing 
the smoothed quantity {S)an, for which semiclassical ap- 
proximations are convergent, one can therefore extract the 
variance of the Si from the scaling in AN of Var[(6>)Aiv]- 
Our starting point for this calculation is 

^4{|^(r)| 2 >A7V = 1± M2kpx)---Im(G osc (r,r,E)) AE 

V 7T 



1 ± J (2k F x) 
Av 



)AE 



(9) 



valid near the boundary, which differs from the expression 
given by Bogomolny [23] only through the inclusion of the 
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FIG. 2: In the semiclassical sums, the four orbits which, as x — > 0, 
coalesce into the same [nearly] periodic orbit, need to be aggregated. 
The top row correspond to two (nearly) periodic orbits. The bottom 
row correspond to two non periodic orbits (p' x = — p x ). 



Bessel function Jo(2k F x) accounting for the Friedel oscil- 
lations. Here, the energy smoothing indicated is similar to 
above except normalized by the energy range, and p° sc is the 
oscillating part of the density of states, given semiclassically 
as a sum over periodic orbits 



p osc (E) 



1 



E 



cos 



'1 2 



|Dct(A/ 7 - 1)| 



1/2 



(10) 



Similarly, the diagonal part of the Green's function 
G R (r, r, E) is expressed as a sum over closed (not necessarily 
periodic) orbits 



G osc (r,r,£;) 



1 1 



ih \j2iixh 



E 



i 



x exp 



-^(r, r, E) — i~T)fj. 



(11) 



In the above expressions, S is the action integral along the 
orbit, r the period, m 12 = 9y final /9pJ / nltlal , M 7 the mon- 
odromy matrix, and <t 7 , rjn are Maslov indices. The tilde on 
the Green's function in Eq. (O furthermore indicates that the 
short orbit giving rise to Friedel oscillations, namely the one 
bouncing off the boundary and returning directly to its initial 
location, is excluded from the semiclassical sum (as in Eq. (0 
it is already taken into account by the Bessel function). 

Inserting Eq. (O into Eq. (0, and performing the integral 
over space, it is important to note that the range of integration 
in the direction perpendicular to the boundary is short, even on 
the quantum scale, and that therefore a stationary phase con- 
dition should be imposed only on the parallel direction. Gath- 
ering however the four orbits shown in Fig |2] and labelling 
their total contribution with the index 7 of the periodic orbit 
to which they converge as x — > 0, one obtains, following the 



usual steps of the derivation of the Gutzwiller trace formula, 
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let 
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where X(6i) is the function introduced in Eq. (O evaluated at 
the angle at which the periodic orbit 7 strikes the boundary 

on the I th bounce, and A 7 ^ f J~f~ l (l ± 



with L ( the total 



v 7 2An^ 

length and n 7 the number of bounces of the orbit 7. 

Eq. ( fT2l can be used directly to compute the Fourier trans- 
form of the Si, yielding a structure very similar to that of the 
density of states Eq. dTOb - Both Fourier transforms can be 
done so as to have peaks with the same shape and positions, 
but different amplitudes. Eq. (fT2l can be furthermore used 
to compute the variance of (S)an- To begin, square Eq. ( [Tot , 
use the diagonal approximation in which orbits are paired only 
with themselves or their time reverse symmetric partner, and 
apply the Hannay Ozorio de Almeida sum rule 112511 



E 

7=p.o. 



Det(Mi-l)| 



= 1 



where only periodic orbits of n bounces from the boundary 
are included. Then for the long orbits relevant here, replace 
the mean length l 1 jn 1 between two successive bounces for 
a specific orbit 7 by d — ttA/C, the mean length between 
bounces averaged for all initial conditions on the boundary. 

Noting (\(9) - A) = 0, this gives 



((S)an-S) = 



2^" 



A(0) - A 



1 

~Kn 



(13) 



As before, angle average (-)g is over the measure d sin 9. 

The structure of this result is quite interesting. The sur- 
viving contribution comes from what may be called the 
"diagonal-diagonal" terms, i.e. pairing not only of the same 
orbit but of the same bounce from the boundary. This 
term contributes only to the variance as it scales as 1/AN; 
see Eq. (0. The terms that would be called "diagonal-off- 
diagonal" give a vanishing contribution. In addition, note the 
absence of a constant term. It implies that the covariance van- 
ishes, which is consistent with our cardioid billiard calcula- 
tions (not shown here). 

Thus, from Eq. ^ and the above considerations 



Var(5 i ) = ^x((A(0)-(A) 9 ) 2 ) 



where 



am - (x) g ) 2 )e 

{ (2 In 2-1) - 
_ 1 (2 In 2 - 1) - 



(14) 
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(15) 
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Comparing this expression with Eq. ©, we see that the semi- 
classical approach and the random plane wave model lead 
to the same result, except for two differences. First, a fac- 
tor two in the prefactor can be traced to dynamical correla- 
tions missed by RPW. The nature of these correlations, which 
are somewhat subtle, will be discussed in 12411 . Second, the 
mean square (X 2 (9))g has been replaced by the variance of 
X(6), giving now a much better agreement with the numeri- 
cally evaluated Si (see Fig.[T}. 

In Eq. (TT2l) . the term proportional to A(0 7 ) can be seen to 
arise from G osc in Eq. (O, while the one proportional to A 7 
originates from p osc . The random plane wave result is thus, 
in some sense, ignoring the latter contribution. As p osc has no 
spatial dependence, its role is not to describe the variations of 
the wave-function, but rather to ensure their normalization (as 
can be seen readily by integrating Eq. (0 over space). There- 
fore the main reason of the failure of the RPW approach is 
due to the lack of individually normalized wavefunctions, and 
more precisely to the fact that Eq. imposes the normaliza- 
tion of the wavefunctions only on average M26J1 . 

Equation ( fT4l i emphasizes the importance of the boundary 
conditions. Indeed, the variance of X(6) is extremely small 
for Dirichlet boundary conditions, which has to be expected 
since the wave-functions are zero near the boundary. As a 
consequence, and as seen in Fig. [T] the fluctuations remain 
smaller than one even for a relatively large number of parti- 
cles, in spite of the linear (k F C) dependence of the variance. 
On the other hand, Neumann boundary conditions yield a log- 
arithmic divergence which can be considered in practice as a 
constant somewhat larger than one. Fluctuations in this case 
are greatly enhanced with respect to the Dirichlet case. 

To conclude, "integrated" wavefunction statistics are intro- 
duced, in part motivated by the need to understand the effect 
of interactions on ground state properties in quantum dots. A 
significant part of their fluctuation properties can be under- 
stood with basic RPW modeling. However, a semiclassical 
framework is developed here, which indicates missing ingre- 
dients of RPW, and in particular that lack of individual state 
normalization leads to a significant overestimate of the fluc- 
tuations. It turns out furthermore that boundary conditions, 
which are often not discussed in the context of mesoscopic 
systems, play an important role. 

We stress finally that only the extreme limit of strongly 
chaotic systems has been treated here. For less developed 
chaos, or systems with some regular dynamics, where some 
eigenstate localization exists, there may be significant en- 
hancements in the fluctuations and correspondingly greater 
effects on ground state properties ll27ll . The consequences for 



experimentally realizable systems and for such systems with 
some form of eigenstate localization is left for future study. 
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